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ABSTRACT 

Cosmological N-body hydrodynamic computations following atomic and molecular chem- 
istry (e", H, H+, H", He, He+, He++, D, D+, H2, HJ, HD, HeH+), gas cooling, star for¬ 
mation and production of heavy elements (C, N, O, Ne, Mg, Si, S, Ca, Fe, etc.) from stars 
covering a range of mass and metallicity are used to explore the origin of several chemical 
abundance patterns and to study both the metal and molecular content during simulated galaxy 
assembly. The resulting trends show a remarkable similarity to up-to-date observations of the 
most metal-poor damped Lyman-a absorbers at redshift z These exhibit a transient na¬ 
ture and represent collapsing gaseous structures captured while cooling is becoming effective 
in lowering the temperature below ^ 10^ K, before they are disrupted by episodes of star 
formation or tidal effects. Our theoretical results agree with the available data for typical el¬ 
emental ratios, such as [C/O], [Si/Fe], [O/Fe], [Si/O], [Fe/H], [O/H] at redshifts z ~ 2 — 7. 
Correlations between HI and H2 abundances show temporal and local variations and large 
spreads as a result of the increasing cosmic star formation activity from z ^ 6 to 2: ^ 3. The 
scatter we find in the abundance ratios is compatible with the observational data and is ex¬ 
plained by simultaneous enrichment by sources from different stellar phases or belonging to 
different stellar populations. Simulated synthetic spectra support the existence of metal-poor 
cold clumps with large optical depth at z ^ 6 that could be potential population III sites at 
low or intermediate redshift. The expected dust content is in line with recent determinations. 
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1 INTRODUCTION 


The current understanding of cosmic structure formation relies 
on gas coolin g and collapse within growing dark-matt er poten¬ 
tial wells te.g. lSaslaw & Zip^ll967l : lGunn & Gotj|l972h . The re¬ 
sulting gaseous fragments represent the sites where star forma¬ 
tion takes place and subsequent stellar evolution processes en¬ 
rich the medium with heavy elements. Metal pollution is deter¬ 
mined by the final fate of stars which synthesise metals during 
their lives and deliver them into the surrounding medium via super¬ 
nova (SN) explosions or stellar winds. The relative production of 
each chemical element is tightly linked to the mass of the progeni¬ 
tor star: massive short-lived stars synthesise and predominantly ex¬ 
pel oxygen a nd the Q-capture elements via type II SNe explosions 
(SN II; e.g. Wooslev & Weaverl 1 19951 : IChieffi & Limongl 1 20041 : 
iHeger & WooslevIbOlOh . whereas less massive long-lived stars can 
enrich the surrounding medium through type la SNe (SN la) which 
produce substantial amounts of iron-peak elements. Therefore, to 
understand the observed abundance patterns at different epochs via 


theoretical models of structure formation it is important to follow 
cosmic chemical evolution of individual elements with appropriate 

delay-times. _ 

Sever al works (see e.g. ICiardi & Ferr^l2005l : lBromm & Yoshidal 
[touI for recent reviews) have shown that primordial structure for¬ 
mation can take place via catastrophic molecular-driven runaway 
collapse in rare high-density peaks already at redshift 2 > 10, 
when the Universe was less than 0.5 Gyr old. At such early times 
metal spreading by primordial sources would rapidly pollute the 
local medium up to metallicitieE Z ~ 10~^Zr^ jMaio et al.l 


I2OIOII2OI lHwise et alj|20ll 

2014; 

Biffi & Maio 

I2OI3) and could 

reach values of ^ 0.1 Zo 

Camnisi et alj|201 1 

; Salvaterra et al.l 


I 2 OI 3 I : IPallottini et alJl2U3 : iMa et al .11 20151) by 2 6, when the 


age of the Universe was roughl y 1 Gyr. At more recent epochs 
larger Z values could be attaine d l lTornatore et al.l2007l : lDave et al.l 
I2OI lL[20T3l : IWuvts et al.ll2014h . This process induced a rapid tran¬ 
sition from a primordial metal-free star formation regime led by 
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^ Metallicities are indicated by Z and are defined as the mass fraction of 
all the heavy elements in the gas phase. 
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non-equilibrium H 2 gas cooling dSaslaw & ZiDoviri967h to an en¬ 
riched o ne characterised by efficient fine-structure metal line cool¬ 
ing (e .g. ISromm & Loebll2003l : ISantoro & Shullll2006l : iMaio et ^ 
l2007h and the subsequent formation of long-lived stars with masses 
< IMq. At lower redshifts, the Universe became increasingly 
dominated by cosmic star formation down to z ~ 2. This 
epoch is responsible for converting a significant amount of gas 
into stars and for a general flattening or decrement of the spe¬ 
cific star formation rate, as observed at redshifts 2 < 2 ; < 6 


(e.g. Daddi et alj2007|; Michatowski et alj2010|; Reddy et al.l2012l : 


IStark et alj2013l : Bouwens et al.l2014l ; Duncan et al.l2014h . Purine 
this period the first intermediate- and low-mass stars evolved and 
locally contributed heavy elements, thereby changing the composi¬ 
tion of the cosmic medium. 

Although the global picture seems qualitatively clear, there are 
still many quantitative uncertainties due to our incomplete knowl¬ 
edge of nucleosynthesis and abundance evolution over cosmologi¬ 
cal times. Additional unknowns are linked to the effects of the envi¬ 
ronment at various z, which might influence the chemical patterns 
of cosmic objects. Despite the utility of simple stellar population 
models in testing typical trends of abundance ratios, they have lim¬ 
ited predictive power in explaining inhomogeneous metal spread¬ 
ing that requires 3D numerical calculations. 

In this respect, a valuable tool for testing chemical abundances 
in the high-redsh ift Universe is offered by damped Ly-a sys¬ 
tems (DLAs) (e.g. lWolfe et alJll986l:lFoltz et al.lll986l : ISmith et al.l 
1986h. both from an observational point of view (see 


Wolfe et al.ll2005aL for a review) and from the theoretical one. Nu¬ 


merical investigations avail able in the literature have mainly fo¬ 
cused on HI propert i es (e.g. Katz_et alJl996l:l Haehnelt et all 19981: 
Pontzen et al.l I 200 A IPopping et al. 20091: Erital et alj 20U : ICenI 


2 OI 2 I: iDuffv et alj|2012l: lYaiima et alj|2012l: iRahmati et alfalSl) 

and neglect heavy-element abundance ratios, which offer a power¬ 
ful tool to study the origin and evolution of cosmic chemical en¬ 
richment. This is one of the primary goals of our paper. 

DLAs are defined as clouds of predominantly neutral gas with 
Hi column densit i es in excess of 2 x 1 0^° atoms cm“^ (e.g. 
iMoller et"^ 1 19981 : IProchaska et al.l l2003h . which appear to be 
associated with galaxies at the faint end of the galaxy lu¬ 
minosity function and represent neutral pr oto-galactic objects 
that are just building up their stellar m ass jFvnbo et al.l l20ld : 
iKrogager et alj|201^ : iBouche et alj|2012h . They are usually cold 
or “warm” structures with typical kinetic temperatu res between a 
few 10^ K jPetitj eanetalj200(j; Srianand et al J201^ and ^ 10^ K 


I Pettini & Cookd 20121 : ICarswell et alj 201^ Noterdaeme et al.l 

2012 h constituted by largely neutra l or molecular hydrogen (e.g. 


Zvizl) constitutea by largely neutra t or molecular nyorogen (e.g. 
Vladilo et alj 1200 ll : ICui et ^ 200^_and with metal content that 
can be measured precisely jPettini ^l2008l:lB^tisti et al.ll2012h . 
Owing to their thermodynamical properties, DLAs are the prefer¬ 
ential sites for the initial stages of gas cooling and star formation 
an d are useful probes of meta l absorption lines at redshift 2 > 
6 (iD’Odorico & Molard l2004l : ISimcoe et al.l l2012l : iBecker et~n] 
I 2 OI 2 I) . These sy stems are also valuable tes t s for stellar evolution 
parameters (e.g. ISalvadori & Ferraral [20l^ : iKulkarni et alj|2013l) 
and chemical evolution mo dels in low-metallicity environments 
or in the early Universe (e.g. Bolton et al.ll20 1 ll : iMaio et '^120131: 
iFinlator et al.ll2013l : lKeating et alj|2014 1. They are precious probes 
of the abundance of nuclides for med during primor dial nucleosyn¬ 
thesis, such as deuterium, D (e.g. lCooke et alj2014h . and have also 
been employed to investigate possible variations in the fundamental 
physical constants (e.g. the ratio between the proton and electron 
mass, /i) at different cosmic epochs (e.g. IWendt & Molardl2012l : 


iKing et ani201 itiRahmani et al.ll2013h . Observationally, DLAs do 
not appear to be prefe rentially associated wi th bright actively star 
forming galaxies (e.g. Fumag^liet^ |2015h. however , both indi¬ 
rect experiments jWolfe & ChenI l200ftjRufeLki_etalJ 201 ih and 
direct searches for host galaxies l lMbller et alj 20041 : Fvnbo et'H] 
l2008h suggest a physical connection between DLAs and star for¬ 
mation. 

Considering the availability of DLA spectra covering a wide range 
of redshifts, we use them in our work as observational tools for the¬ 
oretical predictions of chemical abundance ratios. 

In this paper, we present results for atomic and molecular abun¬ 
dances of cold, gas-rich environments at 2 < 2 < 7 as expected 
from 3D numerical simulations of cosmological structures includ¬ 
ing non-equilibrium chemistry, detailed stellar evolution and chem¬ 
ical enrichment. We compare our theoretical predictions for cos¬ 
mic chemical abundances to the current suite of observational data 
on intermediate- and high -2 observations of metal-poor DLA sys¬ 
tems. We study the effects of cosmic pollution on abundance ratios 
for different heavy elements and investigate the typical trends as a 
function of cosmic time. We find that typical ratios are usually scat¬ 
tered due to simultaneous enrichment by different stellar phases or 
by different stellar populations and are consistent with a general 
decrease of metal content at high 2 . Molecular fractions show tem¬ 
poral and local variations as a result of star formation activity and 
related feedback processes. 

The paper is organized as follows. Details on the numerical im¬ 
plementation and the data sample are given in Sect. chemical 
abundances are discussed in Sect.[3 and the theoretical expecta¬ 
tions from synthetic spectral observations are presented in Sect.|4] 
We summarise and conclude in Sect.O 


2 METHOD 

In the following subsections, we briefly describe the most impor¬ 
tant features of the cosmological calculations we have performed 
(Sect. l2Tt . the selection criteria for our analyses (Sect. lT2t and the 
data samples we have used (Sect. l23] l. 


2.1 Simulations 

The numerical calculations performed in this work extend our pre¬ 
vious studies that investigate the interplay between galaxies and 
the intergalactic medium fr om low redshift (2 < 2) jTescari et al.l 
l2014l : lKatsianis e t ^Jl2015l. A ngus projec t ) to the epoch of reion - 
ization and above 1 Maio et aI]|201(lT201 ll : ISalvaterra et alj|2013h . 
Our hydrodynamical calculations are carried out via the paral- 
lel numerical c ode P-Gadget3, an updated version of P-Gadget2 
(ISpringei2005l) . 


2.1.1 Implementation 

The improved implementation used here combines several physi¬ 
cal processes and, in particular, contains a self-consistent treatment 
of low-temperature cooling by molecules and metals, as described 
in iMaio et al.l ( I2007L 1^1 Oh . We include detailed non-equilibrium 
atomic and molecular chemistry evolution for e~, H, H"^, H~, He , 
He+, He++, H 2 , H+, p, D+, HD and HeH+ jYoshida et al.ll2003l: 

iMaio et al.l2006LE007l : IPetkova & Maicl2012h . determined accord¬ 

ing to the reaction network in Tab. [T] Temporal evolution is com¬ 
puted by solving the following equation for each chemical species 
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Origin of cosmic abundances 3 


Table 1. Reaction network. 


Reactions 

References for the rate coefficients 

H + e- H+ + 2e- 

A97 / Y07 / M07 

H+ + e- -i- H + 7 

A97 / Y07 / M07 

H + 7 -i-H++e- 

A97 / Y07 / M07 

He + e“ —)■ He+ + 2e“ 

A97 / Y07 / M07 

He+ + e“ —r He + 7 

A97 / Y07 / M07 

He + 7 —r He+ + e“ 

A97 / Y07 / M07 

He+ + e“ He++ + 2e“ 

A97 / Y07 / M07 

He++ + e“ —)■ He+ + 7 

A97 / Y07 / M07 

He+ + 7 —)■ He++ + e“ 

A97 / Y07 / M07 

H + e“ —)■ H” + 7 

GP98 / Y07 / M07 

H^ + 7 —>■ H + e” 

A97 / Y07 / M07 

H- + H -r H 2 + e- 

GP98 / Y07 / M07 

H + H+ -).H 2 + +7 

GP98 / Y07 / M07 

H J + 7 ->■ 2 H+ + e- 

A97 / Y07 / M07 

H J + 7 -!■ H + H+ 

A97 / Y07 / M07 

H 2 + + H H 2 + H+ 

A97 / Y07 / M07 

H 2 + H ->■ 3H 

A97 / M07 

H 2 -I-H+ ^H 2 + +H 

S04 / Y07 / M07 

H 2 + e“ ^ 2H + e“ 

ST99 / GB03 / Y07 / M07 

H- +e- ^■H + 2 e- 

A97 / Y07 / M07 

H- + H 2H + e- 

A97 / Y07 / M07 

H- + H+ 2H 

P71 /GP98/Y07/M07 

H- +H+ -7 H 2 + +e- 

SK87 / Y07 / M07 

H 2 + +e- -i-2H 

GP98 / Y07 / M07 

H 2 + +H- ^■H + H 2 

A97 / GP98 / Y07 / M07 

H 2 + 7 —)■ H^ + e” 

A97 / Y07 / M07 

H 2 + 7 -i> 2 H 

A97/R01 / Y03/M07 

D + H 2 ^ HD + H 

WS02 / M07 

D+ + H 2 ^ HD + H+ 

WS02 / M07 

HD + H D + H 2 

SLP98 / M07 

HD + H+ D+ + H 2 

SLP98 / M07 

H+ + D ^ H + D+ 

S02 / M07 

H + D+ -7 H+ + D 

S02 / M07 

D“*" + e~ ^ D + 7 

GP98 

D + 7 ^ + e“ 

GP98 

He + H+ HeH+ + 7 

RD82/ GP98 / M07 

HeH+ + H He + H^ 

KAH79 / GP98 / M07 

HeH+ + 7 -7 He + H+ 

RD82 / GP98 / M07 


Notes: 

KAH79 


1982) 


1997) 


1998); 


7 s tands for ph otons; 
= iKarpas et al ] 


SK87 

GP98 


P71 = jPete^o^^ 

RD82 = iRoberge &J3alg ar^ 


iGalli & Pall j 


19 


19871) ; A97 = 

SLP98 = 
ROl 
S02 


Stancil et al 




Abel et al 


Ricotti et al 


I 


2002 ): 

2003); 


ST99 = lstibbe & Tennys onl j 1 
20011) ; WS02 = IWang & Stancill 1200: 

GB03 = iGlover & Brand j2003h; Y03 = lYoshidaetah _ 

S04 =ISavin et all b004) ; Y07 = lYoshida et alJ i2007h ; M07 = lMaio et alJ 

( 1 ^. 


with number density m at temperature T: 

^ = EE kpq, iiT)n,n,-J2 kii{T)nini, 


( 1 ) 


where t is the time interval, kpq^i{T) is the (temperature- 
dependent) creation coefficient of species i as a result of inter¬ 
actions of species p and q, while kii{T) is the destruction co¬ 
efficient due to interactions of species i and I (see references in 
Tab. [T]). The non-equilibrium treatment imposes additional con¬ 
straints on the simulation timestep, hence abundances are computed 
each ~ 1/10 the hydrodynamical timestep. This choice assures a 
robust conv ergence of our calculations, as previously checked in lit- 
erature (e.g. lAbel et al.ll997Ll2002l ; lYoshida et alj2003l ; lMaio et al] 


1 20071120101 ; IPetkova & Maiol 1201 2l . and references therein), and 
allows us to have reliable estimates for the gas molecular con¬ 
tent even at zero metallicities and in environments where sim¬ 
ple collisional equilibrium is questionable. HI self -shielding cor¬ 
rection s are accounted for on post-processing as in iRahmati et al.l 
) l2013h . taking into a ccount density, temper ature and the evolv¬ 
ing UV backgroun(|3 jHaardt & Madaj[2012l) . We additionally in- 
clude cooling, star formation, stellar evolution and metal spreading 
I Tornatore et alj[200^ ; iMaio et ^120091 l2010l . 1201 ll : iTescari et all 
l2014h by tracing individual elements (He, C, N, O, Ne, Mg, Si, S, 
Ca, Fe) according to available stellar yield s )|Wpqslgy_&(^^ayej 
I 995 I; TOn den Ho ek & Groenewegenl I 1997 I; Wooslev et alJ 2002 


Thielemann et al Jl2Q03ia UV radiation jHaardt & 


ilJ[200l 

tJ(200l[) 


20051; 


and feedback effe cts (see more details in e.g. Springell _ 

ITescari et al.ll2014l) . The UV background can heat the gas and dis 
sociate molecules in the inter-galactic medium (IGM), but cool¬ 
ing and star formation in dense collapsing sites are little affected. 
Feedback from black holes although included i s usually negligi - 
ble at the redshift and masses c onsidered here llBird et al.] 20141) . 
The stellar lifetimes are gi ven bv iMatteucci & Greggid l ll986ll and 
iRenzini & Buzzonil ) ll986l) . The assumed stellar initial mass func¬ 
tion (IMF) covers the ran ge [0.1, 100] Mp and follows the log¬ 
normal shap e proposed by IChabrieil i2003h . numerically approxi¬ 
mated as in iTescari et al.l ( 12014 ). Stars with masses above 8 Mq 
explode as SNe II and inject an energy of 10®^ erg which heats 
the surrounding medium. Lower-mass stars evolve through AGB 
or SNe la phase. Gas and heavy elements can be ejected by star 
forming regions via winds (at 450km/s). We mimic metal diffu¬ 
sion in the interstellar medium by smoothing individual metallic- 
iti es over the neighb ouring particles in the SPH kernel. We refer 
to lMaio et al.l ) l2013ll for extended tests through simulations of iso¬ 
lated objects. 

We adopt a ACDM background cosmological model with present- 
day expansion parameter normalised to 100 km/s/Mpc of h — 
0.704, and baryon, matter and cosmological-constant parameters 
equal to Hp.b = 0.0456, Ho.m = 0.272, Ho ,a = 0.728 
jKomatsu et alj20TIh . Spectral parameters are assumed to be as = 
0.8 and n — 0.96. We consider a box of 12 Mpc//i (comoving) a 
side, where gas and dark-matter fields are sampled by 384^ parti¬ 
cles for each species, resulting in a gas and dark-matter resolution 
of 3.86 X 10® Mq//i and ~ 1.92 x 10® M 0 // 1 , respectively, and 
softening length of 1.5 kpc//i. 

We note that for a precise picture hydrodynamical simulations 
should resolve both rare objects on large scales and very small 
structures on tiny scales. Because of numerical limitations, a fair 
trade-off to assess properly the luminosity functions and halo statis¬ 
tics at high redshift is normally obtained with boxes of ~ 10 — 
30 Mpc//r a side (e.g. iCampisi et al]|201 ll : Isalvaterra et al.|[20l^ ; 
iMaio & Vielll201^ . In terms of space resolution, 1.5 kpc//i soften¬ 
ing length used in this work is enough to pr ovide a realistic de¬ 
scription of galaxies at high and low redshift jTescari et alj|2014l ; 
iKatsianis et alj|201^ . within the limits of the modeling. In terms 
of star formation and metal evolution, our cosmic star formation 
rate densities converge at 2 : < 15 as well as metal distributions and 


^ We note that similar fitting formulas can bear uncertainties due to un¬ 
known ionisation coiTections. 

® We remind the reader that metallicity-dependent iron yields by 


IWooslev & Weaved il995 

) are overproduced by a factor of about 2. See 

e.a.iTimmes et al.Ul995l). 

Francois et alJ j2004|) and references therein. We 


take into account this correction. 
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4 U. Maio & E. Tescari 


filling factors (e.g. iMaio et al.lE)I(il201 iLlillil . These consider- 
ations hold even for different cosmologies and model parameters 
I Maio^Iannu^ 201jJ; Maioet^ 2oI]J^ 201 iL iMaio et aP 
l2012l:lMaio & Khochfarll2012l:lMaio & ViefcOl^ 


It is worth stressing that limitations of stellar evolution results are 
often related to severe unknowns in stellar evolution features. Rec- 
tius, in order to have accurate metal ratios at different epochs it 
is fundamental to follow temporal evolution and mass-dependent 
yields from SN II, AGB and SN la stars. While massive SN II ex¬ 
plosions dominate O and a content in early star forming phases, 
lower-mass AGB stars are particularly important at intermediate 
epochs, because their ejecta contain significant amounts of C, N 
and possibly O. Hence, they can alter the pre-existing chemical pat¬ 
terns. Dwarf stars ending their life as SN la are expected to cause, 
instead, a remarkable boost in Fe content. Since we follow stellar 
evolution from all the aforementioned stellar phases and consis¬ 
tently with mass-dependent yields and lifetimes, our results should 
be very robust. In this work we focus on regular population II stellar 
generations, neglecting the primordial population III (popIII) ones, 
since their contribution to the cosmic star formation is believed to 
be m inor at these epochs (e.g. iMaio et alllioM 1201 ll: Iwise et al.l 
l2014h . 


In addition, when assessing gas molecular content in different en¬ 
vironments it is important to follow detailed molecular evolution 
with a proper non-equilibrium scheme, as the one adopted here. 
Simple phenomenological prescriptions for H 2 and HD can lead to 
misleading conclusions. Effects from dissociating Lyman-Werner 
radiation coming from the first stars can have some impacts on the 
local gas, but the following collapse phases are usually dominated 
by resonant or fine-s tructure transition cooling from newly ejected 
heavy elements (e.g. iMaio et ^l2007h . These have been included 
as detailed in the previous paragraphs. 

Cosmic structures are id entified by means o f friends-of-friends and 
Subfind algorithms (e.e. lDolag et al.ll2009L and references therein) 
and are post-processed to trace masses, positions, velocities, star 
formation rates, H, He, H 2 , HD, C, N, O, Ne, Mg, Si, S, Ca, Fe, etc., 
and all the relevant physical properties in eac h object. To avoid spu¬ 
rious fragmentation or numerical artifa cts ( iBate & Burkerd [19971 : 
iMaio et all l2009l : Ide Souza et al.l l2Q13h . we consider only those 
structures that are resolved with at least 200 gas particles or, equiv¬ 
alently, ~ 10^ total particles. This ensures that behaviour of the gas 
within bound structures is properly treated. 

The cosmological redshifts we consider in this work are 2 = 6.564, 
4.992, 4.186, 2.995, corresponding to cosmic times of ~ 0.85, 1.2, 
1.5, 2.2 Gyr after the Big Bang, respectively (for our adopted cos¬ 
mology). 



Figure 1. Simulated mass function at redshift 2 = 2.995 (p oints) com¬ 
pared to theoretical predictions at the same redshift according to lSheth et alj 
1200 ih (solid line). 



z 


2.1.2 Basic properties 

Numerical results are tested against analytical calculations in 
Fig.m where the mass function extracted from the halo catalogue at 
2 = 2.995 (points) is displayed against corresponding theoretical 
expectations (solid line). The general agreement between numeri¬ 
cal and analytical calculations at all masses is satisfying and sup¬ 
ports the reliability of o ur results. We refer the interested reader to 
our previous works (e.g. lMaio & Viej2015L and references therein) 
for more focused discussions on the effects of different parameters, 
resolution dependence and power spectra in similar simulations. 

In Fig. 12 we also show the evolution of the cosmic star forma¬ 
tion rate density, SFRD, which is co nsistent with observational 
determinations by iDuncan et alj ( |2Q14|) and features an increasing 
trend from ~ 0.02 up to ~ 0.2 M 0 yr“^Mpc“®, as expected at 


Figure 2. Cosmic star formation rate density, SFRD, as expected from 
our work (solid line wit h diamonds) and from observati onal determina¬ 
tions w ith 1 (T error b ars bviReddv & Steidell i2009h (circles). iBouwens et alj 
I 2 OI 4 I 1 (squai'es) and iDuncan et a IJI 2 OIJ) (triangles). 


2^2 — 7, with a flattening around 2 ^ 3 . 

The corresponding metallicity evolution is shown in Fig. [2 where 
maximum, mean stellar and mean gas metallicities are plotted. In 
this figure, Z is computed as the mass fraction derived by C, Ca, O, 
N, Ne, Mg, S, Si, Fe, etc. contributions. The increasing behaviour 
of the mean metal content in cosmic gas is a clear feature of metal 
production from star forming regions. These star forming regions 
can pollute the surrounding gas through feedback effects and can 
boost gas mean metallicities up to 0.1 Zq by 2 < 4. Enrich¬ 
ment can reach maximum values that exceed the solar abundance. 
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0 2 4 6 8 10 12 

X, [Mpc/h] 


0 2 4 6 8 10 12 

K, [Mpc/h] 


0 2 4 6 8 10 12 

X. [Mpc/h] 


Figure 4. Projected physical quantities at various redshift for a slice passing through the center of the box and 0.01 Mpc//i (comoving) thick. The values of 
Xc and j/c refer to the comoving coordinates in Mpc//t along the x and y box sides, respectively. Left, central and right columns refer to mean gas density, 
molecular fraction and metallicity, respectively. 



Figure 3. Metallicity evolution, Z, as a function of redshift, z. Lines re¬ 
fer to mean gas metallicity (solid), mean stellar metallicity (dashed) and 
maximum metallicity (dotted). 


locally. The large spread in metallicity at any z is a consequence of 
two main effects. The first is due to spreading events starting from 


the high-density sites and polluting pristine low-density environ¬ 
ments through wind and SN feedback. The second is due to zero- 
metallicity particles that are particularly abundant at early epochs 
and bring down the mean values when averaging over all the parti¬ 
cles (see further discussion in Sect.|4j. The metallicity of the stellar 
populations (mean stars in the legend) show some temporal varia¬ 
tions, growing from Z ~ 10“^ Zq to Z ~ Zq between z ^ 15 
and z ~ 3, but featuring little evolution already at 2 < 7. 

To have a clear hint on the origin of such metals we will investigate 
in the next sections their abundance patterns, which vary mainly 
with the mass-dependent stellar yields and lifetimes of the actual 
stellar population. 

We warn the reader that, although the overall behaviour should 
remain unchanged, details might depend on the hydro scheme or 
the modeling adopted for feedback processes. They are not ex¬ 
pected to be crucial, though, as previous studies based both on 
SPH and grid codes with different implementa t ions seem to sug¬ 
gest (e.g. Maioet ^ l2010l: iBiffi & Maidl20l3l: IWise et ai]|20l4 


IPallottini et al.ll2014). 


Visual representations of the cosmic objects and filaments are dis- 
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played in Fig.|4] There, we project along the 2 axis density, molec¬ 
ular fraction and total metallicity of a slice passing through the cen¬ 
ter of the box at various cosmic epochs. At 2 = 6.594 primordial 
gas condenses mainly via H 2 cooling which boosts runaway col¬ 
lapse and star formation episodes. Then, massive stars with short 
lifetimes undergo supernova explosions and eject heavy elements 
in the surrounding medium. As a consequence, the Universe can 
be enriched up to metallicities as high as Z ~ 10“® (i.e. ~ 10 
per cent the solar value) already at such epochs. The subsequent 
evolution (e.g. at 2 = 4.992) is characterized by further collapsing 
phases which increases the amount of heavy elements and star for¬ 
mation activity, consistently with Fig.[2 Such events are accompa¬ 
nied by thermal feedback which is responsible for local gas heating 
and partial destruction of existent molecules. This can temporar¬ 
ily halt the SFR as long as H 2 reforms in shocks or in newly col¬ 
lapsing material and metals start cooling efficiently. After a cou¬ 
ple of Gyr, molecular-driven star formation and metal enrichment 
have provided remarkable metal spreading with metallicities reach¬ 
ing values oi Z ^ Zq, as evident from the panels corresponding 
to 2 = 2.995. At all times we recover a sort of typical inside- 
out pattern for cosmic pollution, with heavy elements produced in 
cold dense star forming regions and ejected by feedback effects into 
the lower-density IGM or circum-galactic medium (CGM), which 
is usually characterized by larger entropy and temperature. Envi¬ 
ronmental effects can play a relevant role, since they can dilute 
metallicities (via e.g. shocks by SN explosions or winds) and af¬ 
fect the gas metal content with respect to the stellar one. These 
trends for metal content are consistent with previous stud ies (e.g. 
ISalvaterra et al.ll2013l : iMaio & Vielll2015l : iMa et al.ll2015h . In the 
following we will address these issues in more depth and we will 
investigate the consequences for the resulting elemental ratios. 


2.2 Selection criteria 

The key tool to shed light on stellar evolution and its effects on 
the chemical enrichment of the Universe is the analysis of the 
chemical abundance ratios. In this respect, measurements from 
metal-poor object s at high redshift are available. In particular, 
recent constraints l lCooke et alj|2015l) on the thermal and physical 
properties of metal-poor DLAs suggest typical temperature^ 
of some 10^ K (with a 2a scatter larger than ten thousand K). 
This is consistent with the assumption that DLAs trace a galaxy 
population with low masses and luminositiefl 
Thus, as a first approach in Sect. [3 we will present our theoretical 
predictions by paying attention to those cold/warm gaseous 
structures whose temperature is below 3 x 10“^ K. With this simple 
selection we consider simulated gaseous structures that are, in av¬ 
erage, compatible with the thermal properties of metal-poor DLAs 
and that can reveal typical features of the galaxy formation process 
(Sect. O. The metal content of these systems is characterized by 
[Fe/H] < —2 (see more later). A full discussion about detailed 
distinctions between possible observational classifications (e.g. 
metal-rich DLAs vs metal-poor DLAs) is beyond the goals of 


^ The sample of metal-poor DLAs by ICooke et al.l bOlSl) (their Ta¬ 
ble 4) shows a quite broad temperature distribution with an average of 
960011°°(l<T)iiggo°(2a) K. 

® We note that, although the temperatures mentioned above are compa¬ 
rable to gas temperatures in the Milky Way, the typical properties of the 


2 ~ 3 — 7 gas clouds considered here are significantly different from the 
ones of the Milky Way at 2 = 0 (see Sect.[3). 


Table 2. Reference solar values jAsplund et alj200^ . 


Element: 

H 

C 

N 

0 

Si 

Fe 

Solai' abundance: 

12.00 

8.43 

7.83 

8.69 

7.51 

7.50 


this work. We will rely on metal-poor data simply beca use their 
chem ical content is not affected by dust uncertainties jMolard 
l2006h (for rough estimates see Appendix 1^. Metal abundance 
ratios in these systems are therefore relatively reliable (see next 
Sect. 12.31 . Of course, the one above is a simplified assumption 
that becomes useful when analizing large sets of data and might, 
in some cases, hide systems that could be locally classified as 
absorbing/Lyman-limit systems. 

In order to perform a theoretical analysis akin to observational 
techniques, in Sect. |4] we will generate more sophisticated simu¬ 
lated observations by extracting physical properties along lines of 
sight (LOSs) through the collapsed structures identified with the 
temperature cut. This will ensure that these objects are actually 
DLAs. We will find that the results derived from these two different 
techniques are generally consistent. 


2.3 Observational data 


DLAs are discovered along the line of sight to powerful back¬ 
ground sources (such as quasars, active galactic nuclei, or gamma- 
ray bursts) on the basis of their strong damped Ft I Lya absorption 
line, together with associated metal-absorption lines. 

Since the gas in DLAs is largely self-shielded, most elements 
reside in a single dominant ionisation state in the Hi gas (e.g. 
IVladilo et alj|200lh . all owing the metal content to be derived out 
to redshift 2 ~ 6 (e.g. IPr ochaska et al .11 20071 : IPettini et al.ll200 
lEllison et alj201 lltlBecker et alj2012tl. 

Measurements of molecular fractions are often difficult and com¬ 
plicated by the dissociation ef fects of the ultraviolet background 
in the outer, unshielded zone s jPetitiean et alj[2003 : iLedoux et al.l 
I 2 OO 3 I : iHeinmuller et al.ll2003l . Consequently, mostly upper limits 
are available up to 2 ^ 4 .2 and range between ~ 10“^ and ~ 10“^ 
dNoterdaeme et alj|2008h . Currently, only a dozen measurements 
have been reported, which are typically derived from the coldest 
densest regions of high -2 gas and are generally biased towards 
larger metallicities 1 Noterdaerne_et '^ 120081: ISrianand et af]|20ld: 


lAlbomoz Vasguez et alj2014h. 

In the following, abundance ratios between two arbitrary species X 
and Y are presented with the usual notation: 

[X/Y] = log = log(nx/nY) - logjiix/nv)© (2) 

(nx/riY)© 

where logarithms are base-10, nx and uy are the number fractions 
of the species considered and the subscript © refers to the corre¬ 
sponding solar values. All th e chemical abundance s are scaled ac¬ 
cording to the solar values bv lAsplund et alj ( l2009h . with hydrogen 
normalized to 12 on a logarithmic scale. See Table for the solar 
composition we have adopted in this paper. 

The majority of the data points used in o ur study derive from the 
metal-poor DLA abundanc es reported by Cooke et^ (1201 5h for 
redshifts 2.3 < 2 < 4.5 and iBecker et al.l ( l201lL 2012h forredshifts 
4.7 < 2 < 6.3. First de terminations for high -2 N abundance (a t 
2 = 4.466) date back to lPessauges-Zavadskv et al.l ( l200lLl200% . 
while more recent ones between 2 ~ 2 and 2 ~ 3 can be found 
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Origin of cosmic abundances 1 


in IZafar et alJ i 

2014). Observational data for molecular fractions 

are taken from 

Noterdaeme etalj J2008t) for 1.8 < 2 < 4.2, from 

ISrianand et alj 

2 OIOI 1 for 2 > 3. and fromlAlbomoz Vasauez et al.l 


(l2014h for 2 : ~ 2.66. Further upper limits for H 2 fractions are also 
taken from lNoterdaeme et alj j2008l) for 1.8 < 2 < 4.2. 


We stress that at high metallicities (as large as [Fe/H] > —2) dust 
depletion can affect the gas-phase element abundance ratios, mak¬ 
ing a comp arison betwee n simulations and observations much more 
uncertain ( lMolard[2006h . The advantage of metal-poor DLAs is 
that these concerns are irrelevant and they are a safe and robust 
benchmark to test our numerical predictions. 

The observational samples we adopt throughout this work have 
been collected from a variety of sources and may contain selec¬ 
tion biases. For example, some surveys specifically attempt to find 
the lowest-metallicity DLAs or systems with the highest molecu¬ 
lar fraction. Such biases are very difficult to take into consideration 
and different abundance esti mates can differ by mo re than 1 dex 
(see recent discussions in e.g. iBonifacio et al]|2015lt . For this rea¬ 
son, we will focus mainly on the global trends exhibited by the 
data. Additional constraints on early chemical enrichment might 
be derived by abundance ratios inferred for long gamma-ray burst 
(GRB) host galaxies. Unfortunately, such dat a sample is curr ently 
very poor to draw definitive conclusions (see lMa et alj|201^ and 
will not be considered here. 


3 RESULTS 

In this section we present the main results from our analysis and 
illustrate the evolutionary pathways of cosmic gaseous systems at 
2^2 — 7. Then, we will use measured elemental ratios from metal- 
poor DLAs to probe the chemical enrichment of the Universe in 
addition to the trends of several key abundance ratios over the first 
~ 2 Gyr. 


3.1 Cosmic structures at intermediate redshift 

Fig.|5]presents the basic baryonic properties of our simulated sam¬ 
ple at redshift 2 = 6.594, 4.992, 4.186, 2.995. We show the star 
formation rate (SFR) as a function of gas mass, M, the specific 
SFR as a function of stellar mass, M*, and the stellar fraction, /*, 
as a function of M. 

Typical SFRs are roughly proportional to the gas mass and range 
between ~ 10“® and ~ 10 M 0 /yr. The gas content of our simu¬ 
lated galaxies is up to ~ 10^'’ Mq at redshift 2 ~ 6.594, consisten t 
with previous theoretical expectations at 2 ~ 7 l lMaio et alj|2013h . 
and reaches ~ 10^^ Mq at redshift 2 ~ 2.995, growing by almost 
one order of magnitude in less than 1.5 Gyr. The galaxies that are 
associated with metal-poor DLAs (red asterisks) are normally char¬ 
acterised by low SFR values (~ 10“^ M 0 /yr) with a weak or null 
correlation with gas mass (having usually small Spearman coeffi¬ 
cients, r). We highlight that the estimated S FR for the Milky Way 
is 1 M 0 /yr teobitaille & Whitnevll201(il) . 

The specific SFR, SSFR = SFR/M*, is only slightly decreasing 
with stellar mass, M*, as a consequence of the almost linear de¬ 
pendence of the SFR on and in agreeme nt with observational 
data at redshift 2^2 (e.g. lPaddi et alj|200% . The significance of 
such correlation is quantified by r values that are always around 
~ —0.7 and —0.8. In general, SSFRs are larger at higher red¬ 
shift; at 2 ~ 6.594 there are structures with typical values of about 
~ 10Gyr“^ and in some cases reaching SSFR 10^ Gyr“^. 
This is consistent with estimates inferred from UV emission lines 


at high redshi ft dStark et ^ 201 ^ a nd is in line with r ecent high- 
2 the oretical dSalvat erra et al ]|2()13l ; lDaval et al.ll2013l) and statis¬ 
tical l lde Souza et alj |20l4l20l¥ studies supporting an evolution 
from a very bursty early Universe at 2 > 10 to a more quiescent 
and evolved one at 2 < 6 . At later epochs, when the larger cos¬ 
mic objects have consumed most of their gaseous fuel, the bulk 
of the SSFR distribution drops below a few Gyr“^, with only a 
minor component above ~ 10 Gyr“^, and b roadly agrees with 
the low -2 estimates for 2 ~ 2 — 4 galaxies l lPaddi et al.l 120071 ; 


__ _ 1 P __ 

iMichalowski et alj|2010l ; iReddv et^l2012h . By looking at DLA 
hosts one can see that some of them have large SSFR as a con¬ 
sequence of their extremely poor stellar content, that is often be¬ 
low 0.1 per cent. This is consistent with the simple idea that such 
structures are probably small gas-rich objects that are susceptible 
to bursty epi sodes of star formatio n followed by periods of qui¬ 
escence (e.g. iMarcolini et ^l2006h . As a comparison we mention 
that the Milky Way is instead a gas-poor galaxy with a stellar mass 
of~ 10 “ Mq. 

This is more obvious from the behaviour of the stellar fraction, 
/*, as a function of M. The former is not tightly linked to mass 
(r < 0 . 2 ) and, although larger masses can retain more gas and form 
more stars, feedback mechanisms dominate < 10^ Mq structures. 
They cause a remarkable dispersion due to their more or less severe 
impacts in different environments and affect gas-to-star conversion 
efficiencies in a non-trivial way (see the following subsections). 

B ary on fractions for the Milky Way are difficult to measure. Gas 
content in the solar neighborhood is limited and most of the visible 
mass appears in stars. By assuming a hosting dark halo with a mass 
of~ 10^^ Mq dWatkins et alj 2010 h one can infer a value for /* as 
high as ~ 10 per cent at 2 = 0 . 


3.2 Transient objects 

Besides the general properties of the galaxy population at high and 
intermediate redshifts, it is interesting to focus on some particular 
case studies to investigate the formation of gaseous systems in more 
detail and also to shed light on the evolution properties of DLAs. 
Therefore, we identified a gas clump with a temperature < 10^ K 
and metallicity < 10“^ Zq at redshift 2 ~ 2.995 and traced back 
the constituent gas particles to our higher redshift snapshots. The 
resulting evolution is illustrated in Fig. On the left column we 
plot (comoving) positions at redshift 2 =6.594, 4.992, 4.186, 2.995 
both for the clump particles (big red dots) and for the other particles 
(small green dots) found in the same plane of the center of mass of 
the clump. In the right column we additionally show overdensity 
(5 = p/p) maps with quartile contour levels, as extracted from a 
thin slice containing the reference gas clump. Overall the condens¬ 
ing process of the initially diffuse slightly overdense gas is quite 
clear. It proceeds from scales of some hundreds of comoving kpc 
(i.e. several tens of physical kpc) at 2 ~ 6.594, cools and condenses 
into a denser object whose radius shrinks below a few physical kpc 
at 2 ~ 2.995. The increase of the gas overdensity is well high¬ 
lighted by the color-coded pictures, which show a trend going from 
initial values of 5 1 — 10 up to 5 > 10^ at 2 = 4.186. The corre¬ 

spondence between the red dots on the left panels and the peaks in 
5 on the right panels is striking and consistent with our previous de¬ 
scription. For example, the higher-redshift panels demonstrate only 
moderate density enhancements (5 ~ a few) locally accompanied 
by underdense or void regions in which the lack of material is vis¬ 
ible already from the respective scatter plot. At 2 < 4 the growth 
of the gas clump is in a quite advanced stage: & reaches a few thou¬ 
sands and empty regions are found only in the outskirts (e.g. at 
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Figure 5. The sequence shows the basic baiyonic properties of our galaxy sample at redshift 2 : = 6.594, 4.992, 4.186, 2.995, respectively from top to bottom. 
The three columns refer to star formation rate in M 0 /yr as a function of gas mass in Mq (left), specific star formation rate in Gyr~^ as a function of stellar 
mass in M 0 (center), and stellar fraction as a function of gas mass in M 0 (right). The Speannan correlation coefficient, r, is also quoted on the top-left corner 
of each panel for both the whole galaxy population (empty diamonds) and the metal-poor DLA candidates (asterisks). 
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z= 6.594 





X, [kpc/h] 

z = 2.995 
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Figure 6. Metal-poor DLA formation. The sequence shows the formation 
path of a metal-poor DLA candidate which could be observed at redshift 
z ~ 2.995. On the left column we show the particle distributions of this 
metal-poor DLA candidate at different cosmological times (z = 6.594, 
4.992, 4.186, 2.995, respectively from top to bottom) in comoving coordi¬ 
nates. The right column depicts corresponding overdensity (S) distributions 
(within 10 comoving kpc /h thickness) in 128x 128 pixel maps marked by 
the three quartile contour levels. Color ranges are adjusted to the scale of 
the individual maps to highlight the contrast for a better visualization. 


Figure 7. Metal-poor DLA destruction. The sequence shows the evolution 
and final fate of a metal-poor DLA candidate which could be observed 
at redshift z ~ 6.594. The left column shows the particle distribution of 
this metal-poor DLA candidate at different cosmological times (z = 6.594, 
4.992, 4.186, 2.995, respectively from top to bottom) in comoving coordi¬ 
nates. The right column depicts corresponding overdensity (5) distributions 
(within 10 comoving kpc /h thickness) in 128x 128 pixel maps marked by 
the three quartile contour levels. 


z = 2.995). The entire evolution takes 1.5 Gyr, which is consis¬ 
tent with the cooling timescales expected from the thermal state of 
the gas. The central densest regions would then proceed to collapse 
further to form stars and later suffer from destructive feedback or 
environmental processes. 

This possibility can be verified by following the evolution of a 
pristine DLA which formed at higher redshift. An example is pro¬ 
vided in Fig. [T] where we show the fate of a cold damped system 
forming at z ~ 6.594. This structure hosts star formation at z 5, 
and is significantly affected by stellar feedback processes at z ~ 4. 
At z ~ 2.995 the initial system is almost completely disrupted by 
ongoing SN explosions and thermal heating, which spread the orig¬ 
inal cold gas far beyond the formation site. For this case, the scales 
involved in the process range from some hundreds of comoving 
kpc (i.e. tens of physical kpc) at z ~ 2.995 down to a few physical 
kpc at z ~ 6.594. The existence of such dense objects {S > 10^) 


at these redshifts is a clear hint of primordial structure formation 
led by pristine molecular formation in the first Gyr. Star formation 
feedback (e.g. SN heating, shocks), partly accompanied by envi¬ 
ronmental effects (gas stripping, interactions) at later times, plays a 
major role in shaping the gaseous systems. 

Evidence in support of this picture is provided by the panels in the 
second row (z = 4.992), where typical <5 values slightly decrease 
and gaseous material starts being ejected due to the first star forma¬ 
tion episodes. At later times, the amount of gas that is expelled is 
higher (see third row) and the density in the central region of the 
clump drops down to 5 ~ 10^ — 10^. Simultaneously, structure for¬ 
mation in nearby sites sets in and further contributes disrupting the 
system by tidal interactions. In the bottom panels, filamentary-like 
structures are recognizable, but the constituents of the initial clump 
(big red dots) are almost completely scattered around. 

These two examples suggest that the cold/warm structures seen as 
metal-poor DLAs can possibly be transient cosmological objects. 
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These would represent the final stages of gas collapse and frag¬ 
mentation, captured when gas cooling is bringing temperatures to 
low values and locally the gas is condensing to high overdensities 
(5 > 10^). Their inner regions will eventually form stars and the 
associated feedback mechanisms will evacuate partially, if not com¬ 
pletely, the system’s gas reservoir. 

A more quantitative statement can be done by checking for the 
cooling capabilities of cosmic gas. The cooling time, tcooi, is gen¬ 
erally given by 

E 

tcool = -:- (3) 

E 


Here we will largely focus on the theoretical gaseous chemi¬ 
cal patterns at such epochs. Previous results have shown that 
typical Z values in the first Gyrs are distributed in a range of 
logic (Z/Zq) between —1 and —3.5. Rare outliers at lower Z can 
be present although with probabilities below 0.1 p er cent. We re¬ 
fer the interested reader to our previous works (e.g. Camgisie^^ 
201 lUSalvaterra et alj|2013l: lOaval et al.ll2013l: iMaio & Vielll2015l: 


Ma et alj2015t) for more detailed analyses about probability distri¬ 


bution functions of metallicities at different redshifts. 


3.3.1 Molecular and metal content 


where E and E are the thermal energy and the energy loss rate. 
This equation can be simplified into 


_3 kgT 
^ 2 A riH 


(4) 


where the prefactor 3/2 holds strictly for mono-atomic gas, ks 
is the Boltzmann constant, T is temperature, A is the metal- and 
temperature-dependent cooling function, uh is the H number den¬ 
sity. For typical values in cgs units we can write 


fcooi ~ 2 X 10 


10 (T/IO") 


(A/10- 


riH 


10'" 


(T/10^ 


(A/10-22) nn 


yr (5) 


and easily get some useful estimates. Given that the cooling times 
for diffuse gas (~ 10“® cm”®) are of order ~ 1 Gyr, much longer 
than catastrophic molecular-driven runaway cooling (<C 1 Gyr), it 
is most likely to observe a metal-poor DLA prior to the onset of star 
formation (Fig. [6j, before the effects of stellar feedback can heat or 
evacuate the hosted material (Fig.|7](. 


3.3 Chemical patterns 

The most reliable probe of cosmic chemical evolution is the study 
of abundance ratios for different chemical elements. Abundance ra¬ 
tios depend primarily on the stellar yields, timescales and IMI0, 
whereas the absolute gas-phase metallicity depends on the condi¬ 
tions of the environment, the star formation history, as well as feed¬ 
back mechanisms. We will test chemical evolution during the first 
Gyrs by means of molecular fractions and several important abun¬ 
dance ratios for the most common heavy elements. In particular, 
among the tracked species, we will consider H 2 molecules and the 
main species produced by high- and low-mass stars, such as C, N, 
O, Si and Fe, for which reliable observational data exist. Our refer¬ 
ence indicators are [O/H] and [Fe/H], which are extremely useful 
to trace different stellar progenitors, since O (and the a elements) 
are mainly produced by short-lived massive stars, while Fe is more 
abundantly produced at later times by long-lived low-mass stars. 
Our findings are summarised in Figs. ElEBQiil and where we 
plot the gas mass, M, average molecular fraction, Xmoi, metallicity 
over H, Z/H, and a number of element ratios for all of the simu¬ 
lated galaxies (empty diamond). We highlight the metal-poor DLA 
host candidates with red asterisks. As a comparison, we overplot 
observational measurements with error bars or experimental limits, 
as described in Sect. 12.31 We note that, since we are at relatively 
high or moderate redshifts (z > 3), stellar evolution and metal en¬ 
richment are still in their infancy (the corresponding cosmic time is 
<2 Gyr). 


® Judies bv IConrov & van DokkumI l2012h (but see also lLa Barbera et ^ 
I 2 OI 5 I for different interpretations) suggest further IMF dependencies on 
[Mg/Fe], which is a proxy for a elements. 


The first two panels of Figs. 1811 D show the metallicity and molecu¬ 
lar content in our simulated galaxies for different gas masses (as in 
Fig.|5}, while the third panel displays the Xmoi — [Z/H] trend. 
Mean molecular abundances at ^ = 6.594 for the simulated 
galaxies span roughly two orders of magnitude, ranging between 
~ IQ-S-S ^ 10“®'®, while metallicities are between [Z/H] ~ 
—3.5 and ]Z/H] ~ — 1 (Fig.[^. At later times the molecular con¬ 
tent ranges between 10“® < Xmoi ^ 10“® and the metal content 
between —4 < [Z/H] < —0.5 (Fig.[TTJ. The behaviour of chemi¬ 
cal content with gas mass highlights the capability of larger objects 
to retain more metals and host molecule reformation thanks to the 
deeper matter potential wells. Lower-mass structures, instead, are 
more strongly affected by star formation and feedback events and 
their gas gets easily heated by nearby star forming regions. Conse¬ 
quently, molecules are dissociated and the resulting fractions cover 
a broad range around Xmoi ~ 10“® — 10““* (Fig. i9Tni for gas 
masses below ~ 10® M©. A similar evolution in low-mass objects 
is recognizable for [Z/H] as well, while in the massive end [Z/H] 
tends to values of ~ —1 (Fig. I8I9I ) or even higher (Fig. fToTTB . 
High metallicities are additionally found in dense clustered regions, 
quite independently from the gas mass, where star formation takes 
place efficiently (thanks to the locally high molecular content) and 
spreads metals in the sun'oundings through an inside-out mode. 

As a consequence, objects with large Xmoi and high metal content 
(upper region in the Xmoi — [Z/H] panel) are those that are hosting 
star formation and have already experienced heavy-element pol¬ 
lution, while objects with large Xmoi and small Z values (lower 
region in the Xmoi — [Z/H] panel) are those that are potentially 
star forming, but the onset of star formation has been quite recent, 
hence metal spreading has not yet been able to enrich the medium 
significantly. We note that Xmoi — 10“® is roughly the amount 
of primordial H 2 that forms soon after the decoupling epoch (i.e. 
an approximate indicator of primordial Xmoi), therefore higher val¬ 
ues actually testify the collapse of gas in the dark-matter potential 
wells that boosts H 2 and HD production. On the other hand, values 
lower than 10“® reflect environmental effects that determined 
molecule dissociation around ~ 10^ K (as clearly visible at lower 
redshifts). 

These considerations are useful to understand the features of 
molecular and metal content at later times, when star formation 
and feedback effects become progressively more important. They 
are responsible for increasing the scatter of Xmoi and Z values be¬ 
tween 2 = 4.992 and z = 2.995 by several orders of magnitude, 
due to ongoing stellar evolution and feedback mechanisms that con¬ 
tinue to enrich the surrounding medium and spread heavy elements 
over increasingly larger volumes. As a result of the mechanical and 
thermal feedback from SN explosions, molecular fractions drop 
down to ~ 10“®. From the Xmoi — [Z/H] panels in Figs. [8}{^ 
it is clear that DLAs with low molecular content and low metallic¬ 
ities are usually severely affected. Such structures are likely to be 
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Figure 8. Chemical abundances for the entire simulated population (black diamonds) and that hosting metal-poor DLAs (red asterisks) at redshift 2 ~ 6.594. 
The first row displays molecular (first panel) and metal content (second pa nel) as a function of gas mass, metallicity vs molecular fraction (third panel) and 
[Si/0] vs. [C/0] ratios with data points and eiTor bai's (fourth panel) from iBecker et alj feOlA . The second row collects important element ratios, such as 
[0/F e] vs. [Fe/H] ( first p anel), [Si/Fe] vs. [Fe/H] (second panel), [C/0] vs. [O/H] (third panel), [N/O] vs. [O/H] (fourth panel), with observational data taken 
from iBecker et alJ i l2012h (see text). Contour le vels refer to ia and 3a confidence. Shaded areas in the bottom row refer to the range allowed by observational 
constraints. Reference solar data are taken from lAsplund et al] i2Q09l) . 


2 = 4.992 




[C/0] 



Figure 9. Same as Fig.O Simulation snapshot refers to z 


4.992. [N/O] data point (light blue) is from lDessauges-Zavadskv ^ l^ i200i|) . 


externally polluted by nearby galaxies and suffer molecular disso¬ 
ciation from the UV background (compare the third panel of Fig.[^ 
and the one of Fig. [ID- Available data points at 2 = 4.186 and 
at 2 = 2.995 show agreement for the Xmoi — [Z/H] distribution, 
with the caveat that observed values at Xmoi ~ 10“® — 10“^ refer 
to data that probe the cold neutral medium of high-redshift galax¬ 
ies (with temperatures ~ 80 K) and do not refer to the average 
content of the whole hosting structures, for which the molecular 
fraction would be significantly lower. The range of values are nev¬ 
ertheless in broad agreement. In fact, most of the DLA gas is diffuse 


and warm (T > 3000 K), consequently, H 2 is scarce and detected 
only in ~ 10 per cent of the cases (e.g. lAlbomoz Vasguez et al.l 
l2014h . Additi onal upper limits for Xm oi with available metallic¬ 
ity estimates (iNoterdaeme et al.ll2008h shown in Figs. llOIil II lie 
at Xmoi < 10“"^ and are in agreement with the simulations. We 
note t hat high -2 observational sam ples (e.g. iRafelski et alJl20iT 
|2014 [Jorgenson et al.ll20lT l2014ll usually consist of objects with 
[Z/H] > —3. This floor is quite interesting, because when com¬ 
pared to stellar metallicities it is found that values inferred from in¬ 
dividual stars do not feature such a plateau and attain values lower 
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z = 4.186 
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Figur e 10. Same as Fig.[^ except the simul ation snapshot is taken at redshift 2 : ~ 4.186. Data points for the metal ratios (blue) are taken fromICooke et al, 


20^ ). [N/0] data point (light blue) is from[Dessauges-Zavadskv et alj j200l[). Data points for molecular fr actions (cyan) are taken from lNoterdaeme et al, 
200^ ) (1.8 ^ z < 4.2). ISrianand et aTIfeOlol) lO"*^ al z > 3) and rAlbomoz Vasquez et anj201‘^ — 10~^'^ at z 2.66) and refer to the 


coldest molecular components within the hosting DLA. Upper limits for H 2 fractions (green) are taken from lNoterdaeme et al.l i2008h (1.8 ^ z < 4.2) 


2 = 2.995 



Figur e 11. Same as Fig.[^ except the simulation snapshot is taken at redshift 2 : ~ 2.995. Data points for the metal ratios (blue) are taken from ICooke et alj 
i2015l) . while data points for [N/0] ratios (light blue) are from IZafar et al.l j2014h . Data points for molecular fractions are the same as previous Fig. Eo] 


than —4 (e.g. lCavrel et ^l2004l : ICaffau et al.l201ll : iRoederer et al.l 
|2014 . Considering that it is quite difficult to measure signals from 
low-metallicity gas, the Xmoi — [Z/H] panels suggest a non-zero, 
albeit small (below ^0.1 per cent), probability of finding them as 
a result of the metal enrichment process from star forming regions 
towards the outskirts of galaxies and beyond. 

3.3.2 Abundance ratios 

We continue our discussion on metal enrichment from stellar evo¬ 
lution by showing selected relevant abundance ratios as functions 


of both [Fe/H] and [O/H], In fact, [Fe/H] and [O/H] can be used to 
track the temporal evolution from early times, when massive SNII 
progenitors eject a lot of O, to later times, when long-lived low- 
mass stars eject large amounts of Fe. This means that going from 
low to high [Fe/H] is like following stellar evolution from early to 
later stages or, equivalently, from short-lived larger stellar masses 
to long-lived smaller stellar masses. We cross-check our results by 
looking at [0/Fe] vs. [Fe/H], [Si/Fe] vs. [Fe/H], [C/O] vs. [O/H], 
[N/O] vs. [O/H] and [Si/0] versus [C/O] (see panels in Figs. lSinT l. 
For the various distributions we quote the mean and Icr dispersion 
in the corresponding legends. We now describe our results on gas 
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metallicitie|3 and discuss more in depth their physical origin. 

In general, there is some redshift evolution in the various ratios 
with variable spreads due to contaminations from chemical feed¬ 
back by different stellar phases. This is quite remarkable, because 
the scatter in the ratios highlight the limits of the closed-box model 
and possible deviations from a single stellar population scenario or 
simple semi-analytical expectations. 

More in detail, at early regimes C, O and a elements are abundantly 
produced by massive stars and their trends can be constrained by 
higher-redshift {z > 5) data on [Si/O] versus [C/O]. Values of 
[Si/O] are nearly solar already at 2 = 6.594 and lie with in the 
error bars suggested by the observations teecker et al.ll2012ll . This 
is not very surprising, because a elements are produced in SN II 
by the same nuclear process (a capture) starting from O nuclides, 
hence their abundance is usually proportional to the one of O. As 
a result, mean [Si/O] values for the whole (metal-poor DLA) pop¬ 
ulation evolve only from —0.07 (—0.08) at 2 = 6.594 to —0.03 
(—0.05) at 2 = 2.995 reaching —0.06 (—0.06) at 2 = 4.992 and 
—0.05 (—0.08) at 2 = 4.186. The mean values for [C/O] instead 
increase from ~ —0.4 up to ~ —0.3, being more sensitive to the 
ongoing stellar phases - see also the next discussion. 

We then check the behaviour of [O/Fe] as a function of [Fe/H], 
which represents a very important diagnostic of SN II and SN la 
evolution. [O/Fe] predictions lie in a range around mean values of 
~ 0.4 at 2 = 6.594 and 2 = 4.992 and they are led by short-lived 
massive SNe that eject large amounts of O. At lower redshift a knee 
around [Fe/H]~ —2 is recognizable and marks the transition from 
short-lived massive SNe II to long-lived low-mass SNe la that pro¬ 
duce huge amounts of iron and cause a decline of [O/Fe] values at 
2 = 4. 186 and 2 = 2.995 . Si mulated data are in a greement with 
data bv ISecker et al.l l l2012h and ICooke et alj l l2015l) at 2 = 6.594 
and 2 = 4.992 (despite the lack of observational information on 
[Fe/H]) and with ICooke et al.l ( 1201 Sh at 2 = 4.186 and 2 = 2.995. 
Also in this case we note a large spread of the simulated data which 
is a consequence of the mixing of different stellar phases due to on¬ 
going star formation and feedback mechanisms. 

Similar considerations can be drawn about [Si/Fe] ratios, that 
demonstrate super-solar abundances and are driven by O evolution 
during their early phases and by Fe production at later times. 

The [C/O] ratios for the whole galaxy (metal-poor DLA) popu¬ 
lation present typical mean values that are slightly sub-solar, be¬ 
tween —0.41 (—0.33) at 2 = 6.594 and —0 .31 (—0.28) a t 2 = 
2.995 (consi stent with e.g. lPettini et al.|[200l : lFabbian et al.ll2009l : 
ISecker et al.ll2012l) and feature a U-shaped trend with a minimum 
of ~ —0.5. Hints of this trend can be already seen at 2 = 6.594 and 
2 = 4.992, when primordial pollution episodes by massive SNe 
eject the first metals and pollute the medium. At later times, [C/O] 
can reach the solar value (i.e. at 2 = 4.186) or even overcome it 
(2 = 2.995). In general, the relative increase of carbon abundance 
[C/O] is due to the metallicity-dependent carbon yields of massive 
stars and t he delayed release of C from low- and interm ediate-mass 
stars (e.g. lAkerrnan et al.l l2004l : ICescutti et al.l 1 20091) . Therefore, 
due to early SN enrichmen0[O/H] increases up to about —2 or 
— 1 while [C/O] decreases down to ~ —0.5 (Fig.[8]and|9](. After 
a few 10® yr, AGB stars expel further heavy elements and make 


^ In the following analysis we will refer to gas metallicities, but it should 
be noted that the majority of metal poor stars are substantially o-enhanced 
and metallicities measured via a elements in DLAs and via Fe in stars are 
not totally consistent. 

® Pristine or low-Z SN II have yields that account for values of [C/O] 
typically larger than the ones derived by high-metallicity SN II. 


the [C/O] trend at [0/H] > —2 increase (Fig.fTOlandfTTT). At low 
metallicities ongoing explosions of massive SNe determine the in¬ 
crease of [C/O] ratios. Observations agree with theoretical expec¬ 
tations quite independently from 2 or from unknown [O/H] values. 
Uncertainties appear for N abundances, both theoretically and ob- 
servationally. Theoretical limitations are not due to major failures 
in the calculations, but to the specific input yield tables. N evo¬ 
lution is in fact a very problematic issue for which there is still a 
lack of understanding and a missing agreement about its production 
processes. These are due to non-trivial dependencies on several pa¬ 
rameters of nucleosynthesis and stellar-structure models, such as 
met allicity, magnetic fie l ds, rotation, winds, binary systems, etc. 
(seelFra nqois et ahl 20041: Pettmiet^D 200^ Cescutti & Chiappinil 


l20ld : Cescutti et al.l 20131 : Zafar et ^ 2014 . for further details). 


Furthermore, most of the available data give only upper limits to 
[N/0] that qualitatively point towards an increasing (almost bi- 
modal) spread at low metallicity. Overall, the broad [N/O] values 
suggest that metal pollution derives from stars forming in quite dif¬ 
ferent metallicity environments, from very metal-poor (as testified 
by low [N/O] ratios) up to very metal-rich ones (as in the cases 
with high [N/O]). Considering that N is poorly produced by mas¬ 
sive stars and is mostly expelled by intermediate-mass stars during 
their AGE0 phases it is possible that when intermediate- or low- 
mass star formation took place (at 2 < 7) the cosmic medium could 
be locally enriched up to solar levels with a spread in the ratios of 
a few orders of magnitude. Such large dispersion in [N/O] ranges 
from solar values down to [N/O] ~ —1.5 and is consistent with 
the delayed deliv ery of primary N into the surround ing medium 
with respect to O dPettini et alj2008l : IZafar et alj2014h . In addition 
to the few data points available for high -2 meal-poor DLAs, several 


determinations of extra-galactic H II regions (e.g. 

Izotov & Thuad 

I 2 OO 4 I; Ivan Zee & HavneJl2006l: IPettini et al.ll200S 

, and references 


therein) give [N/O] ratios between —2 and ~ 0 and feature a de¬ 
creasing trend from values in the interval ~ [—1,0] at [O/H] 0 

down to [N/O] ~ —2 at [O/H] ~ —1.5. At very low metallicities 
observations seem to suggest a sort of plateau for [N/O] (albeit very 
scattered), du e to the primary production of N that basically traces 
O abundance dPettini et alJbOOsI : IZafar et alj|2014h . Given the ob¬ 
servational difficulties and the current theoretical limitations, more 
studies and improved yield calculations are needed in this respect 
in order to reach a general consensus. 

In conclusion, we notice that theoretical ratios are generally within 
the observational range and, despite the relatively large error bars 
of the available data, this is extremely encouraging. 


3.3.3 Cosmic chemical evolution 

We now explore the cosmological evolution of metal enrichment by 
comparing the elemental ratios of the DLA population at different 
redshifts. We note that, at the epochs considered here (2 ~ 2 — 7), 
the Universe is already a-few-Gyr old and hosts numerous SN II ex¬ 
plosions and AGB pollution events. Long-lived SN la are still rare 
at 2 ~ 5 — 6, but become progressively important at later times. 
To address the implications of different stellar evolutionary phases 
we consider some indicative elements, such as C and O for SN II 
and AGB, Fe for low-mass SN la and Si as a proxy for a-element 
production by massive stars. 

In Fig. [T^ we collect the expected values of [C/O] (upper panel). 


^ Model unc ertainties for AGB yields aie still large, though. See e.g. 
iKarakaj l201(t) for a comprehensive discussion. 
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redshift 



redshift 



redshift 

Figure 12. Observational determinations (blue crosses) of abundance ra¬ 
tios at different redshifts for [C/0] (upper panel), [0/Fe] (central panel) and 
[Si/Fe] (lower panel) with la error bars compared to our theoretical predic¬ 
tions (red solid line with asterisks) with corresponding 1 a dispersion. 


[O/Fe] (central panel), [Si/Fe] (lower panel) with Icr error bars 
for redshifts z > 2 and compare these to the observational data 
points. In general, our numerical calculations agree with obser¬ 
vations, since the theoretical predictions are within the error bars 
or consistent with the experimental limits for all redshifts. Values 
referring to the low-temperature population and the whole galaxy 
population at different redshifts are quite similar (further details are 
in Sect.[333J. 

Theoretical [C/O] ratios take values of —0.28, —0.31, —0.32 and 
—0.33 for the metal-poor DLA population and —0.31, —0.36, 
—0.39 and —0.41 for the whole galaxy population, at red¬ 
shift 2.995, 4.186, 4.992, 6.594, respectively. These estimates 
are in agree ment with SN II yields for different initial stellar 
metallicities llWoo slev & Weaved 1 19951 : IChieffi & Limongil 1 20041 : 
iHeger & Wooslev 2010h and with most of the data. Interestingly, 
we see a spread in the simulated data that spans a factor of a 
few. For the metal-poor DLA sample, the theoretical dispersions 
are ~ 0.20 at 2 ~ 6.594 and 0.19 dex at 2 ~ 2.995, over 
a timescale of ~ 1.5 Gyr. However, values of 0.13 and 0.12 are 
found at 2 = 4.186 and 2 = 4.992. This spread is a sign of 


enrichment by multiple star formation episodes which take place 
in different environments and which pollute the medium through 
stars in different evolutionary stages and with different metallici¬ 
ties. Indeed, the ratio between C and O can easily deviate from the 
single-stellar-population expectations for a number of reasons. For 
example, enrichment by SN II belonging to the same simple stellar 
population (i.e. born in a medium with similar conditions, with sim¬ 
ilar initial metallicities and at similar epochs) would feature a given 
[C/O] at any given time. However, if SN progenitors were bom si¬ 
multaneously in media with different initial metallicities then the 
SN II yields would differ as well and the resulting abundance ra¬ 
tios would change accordinglvF^. In the same spirit, if, c. p., stars 
were not bom simultaneously the various stellar evolution phases 
would proceed in an asynchronous way and SN II explosions could 
pollute the surroundings when AGB stars are also spreading heavy 
elements. Consequently, abundance ratios would deviate from the 
predictions for a simple stellar population, because of the presence 
of multiple populations. Such deviations would be more important 
mostly at later times, when the evolved stars that coexist with short¬ 
lived stars become effective at polluting the surrounding medium. 
In fact, although at 2 ~ 2 — 4 SN II are still dominanj^, SN la and 
AGB stars can contribute already at 2 < 6, when the cosmic time 
elapsed is around 1 Gyr or larger. 

Obviously, these considerations apply to all the elemental ratios and 
to [O/Fe] in particular. Given the larger error bars, the [O/Fe] evolu¬ 
tion is essentially flat with typical values around 0.4 at all redshifts. 
The galaxy population features [O/Fe] values of 0.35, 0.40, 0.42 
and 0.45 at 2 = 2.995, 4.186, 4.992, and 6.594, respectively, while 
the metal-poor DLA candidates have corresponding ratios of 0.32, 
0.42, 0.38 and 0.41. Such results are consisten t with recent obser¬ 
vational determina tions of metal-poor DLAs dCooke et alj|201 j : 
iBecker et al.l201 j) and the error bars are consistent with the whole 
data sample at any 2 . The slight decrement in [O/Fe] at 2 ~ 3 — 4, 
when the cosmic time is t{z) ~ 2.1 — 1.5 Gyr, reflects the ad¬ 
ditional Fe contributions from low-mass stars that formed roughly 
~ 1 Gyr earlier, at 2 6 — 10, and that explode as SN la enrich¬ 

ing the medium at these very epochs. The interplay with coexisting 
short-lived stars, just forming in recently collapsed regions, induces 
the large dispersion in the abundance ratio at 2 ~ 3. The trend in 
the expected spread is not regular, but rather stochastic and ranges 
within a factor of ~ 4 so that the standard deviation, a, results to 
be about 0.15 dex at 2 ~ 4.992 and 0.65 dex at 2 ~ 4.186. 

Finally, the whole galaxy population at 2 = 2.995, 4.186, 4.992, 
and 6.594, is characterized by [Si/Fe] values of 0.32, 0.35, 0.36, 
0.38, respectively. There is very little difference from the metal- 
poor DLA population, whose corresponding [Si/Fe] ratios are 0.27, 
0.34, 0.32, 0.33. Similarly to [O/Fe] also [Si/Fe] features r oughly 
constant valuej^ of 0.3 (consistently with e.g. Iviadild[2003 : 
IWolfe et al.ll2Q05a) with a decrease at 2 ~ 3 — 4 and a larger scat- 


Large deviations would be expected for several elemental ratios, such 
as [N/0], [Al/O], [Fe/0], etc.. 

We remind the reader that at these epochs the cosmic star formation 
rate density reaches its peak values. 

Here we should keep in mind that a given metal-poor DLA might have 
had a different star formation history than a typical DLA. Similarly, a typ¬ 
ical DLA may have had a different chemical evolution to the Milky Way. 
Based on studies of the local-group dwarf galaxies, the Milky Way and 
the dwarf spheroidals seem to deviate from this plateau in [a/Fe] when 
[Fe/H] ~ —2 roughly. Thus, it is interesting that the red contours in the 
ratio distributions (Figs. [sTTB tend to move to lower values than the black 
points with decreasing redshift, which could be exactly this effect. 
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Figure 13. Mean abundance evolution for the whole simulated [C/0] sam¬ 
ple (black solid line with diamonds) and for the metal-poor DLA candidates 
(red solid line with ast erisks) compared to available observational data by 
ICooke et"^ i20 1 Stl and iBecker et alJ i2012h . The blue dotted line marks the 
observational average over the whole redshift range with corresponding 1 a 
dispersion (green shaded area). The blue dashed line marks the observa¬ 
tional average refemng to the redshift range 2.9 < 2 < 6.7, that is directly 
comparable to our theoretical expectations, with corresponding Icr disper¬ 
sion marked by the blue shaded area. 


DLA sample (red solid line and asterisks). The theoretical expecta¬ 
tions for the whole galaxy population are typically below the ones 
referring to metal-poor DLAs, because these latter sample mainly 
the low-metallicity end, hence the resulting average ratios are a bit 
larger. Slightly lower mean values are obtained when restricting the 
observation range from z > 2 to 2.9 < z < 6.7. This is due to the 
[C/O] values observed at z ~ 2.4 that increase the mean of ~ 0.05 
dex, i.e. about 12 per cent. We note that all the data are compatible 
with metal spreading from regular population II stellar generations. 
The effects of SN la going off at z < 4 — 5 should result in larger 
line widths as a consequence of the establishment of a more and 
more caothic environment. 


4 OBSERVATIONAL SIGNATURES 


As a final task, we compute simulated spectra (Sect. I4.lt in or¬ 
der to retrieve further information that are directly comparable 
against observations. We generate simulated observations by ex¬ 
tracting spectra of Hi, H 2 and metals along lines of sight (LOSs) 
th rough the cosmologi c al box and fol l owing the pr ocedure outlined 
in iTheuns et al.l ( 1 199^ : Tescari_et_alJJ200^ 201 ih . We model self- 

alj 


shielding by following iRahmati et alj i 2013 ) (their fitting formula 
is valid up to z ~ 5, even though we assumed its validity up to 
z = 6.594). We then discuss the column density distribution func¬ 
tions (Sect. l4^ and the statistical trends derived for different red- 
shifts (Sect. l4.3l l. 


ter at z ~ 4.186. In fact, standard deviations vary between ~ 0.10 
(z = 4.992) and 0.46 (z = 4.186). Also in this case, as for 
[O/Fe], the drop at z ~ 3 — 4 is the result of SN la enrichment that 
boosts Fe content with respect to O and a elements. 

We note that data in Fig. [I 2 ] can be easily explained by metal 
spreading and feedback effects from regular stellar populations at 
“2 z <7 without invoking exotic popIII generations. This is con¬ 
sistent with several studies in the literature predicting a rather small 
popIII contribution to cosmic star formation. 

We finish by stressing that evolution and abundances of different 
heavy elements are strictly linked to timescales and yields of dif¬ 
ferent stellar phases, which also depend on the patchy distribution 
of metals in the original star forming sites and on the effects of 
feedback mechanisms (such as inflows and outflows). This poses 
serious limits on homogeneous models or on the closed-box model 
and highlights the need of detailed 3D modelling, as we reported 
herein, to predict both the trend and the dispersion that are observed 
for elemental ratios at various redshifts. 


3.3.4 Direct comparisons 

In the following, we briefly check the implications of outliers in 
a more direct comparison between observations and theory for the 
case of [C/O] ratios. Fig.[T^shows mean abundance evolution for 
the whole simulated [C/O] sample (solid black line) and for the 
DLA candidates (solid red line) compared to available observa¬ 
tional data. The horizontal dotted line and the light-shaded area 
mark the mean value of the whole observational data sample and 
corresponding la deviation. The horizontal blue dashed line and 
the blue-shaded area mark mean value and la deviation for data 
in the redshift range 2.9 < z < 6.7. This latter is more directly 
comparable to the theoretical expectations for our galaxy popula¬ 
tion (black solid line and diamond symbols) and for the metal-poor 


4.1 Spectral analysis 

Fig. [T 4 I shows some physical quantities interpolated along three 
LOSs around the most massive halo at z = 6.594. The black solid 
line is a LOS through the center of mass (CM) of the halo. The red 
dashed and blue dotted lines are two LOSs at, respectively, 50 co¬ 
moving kpc/h and 100 comoving kpc//i from the CM in the same 
direction. 

In the top panel the density of the gas (in nn/cm^) is plotted. A 
sharp spike corresponding to the position of the center of the halo 
is clearly visible at u « 1000 km/s. Apart from the central spike, 
the gas density profile does not change remarkably when consider¬ 
ing LOSs far from the CM. 

The same is true for the temperature profile (second panel). We 
note that the temperature profile is characterized both by hot ma¬ 
terial with temperatures above ~ 10"* K and by cold material with 
values below ~ 10^ K. The former one represents turbulent gas 
that has been heated mainly by stellar evolution and star formation 
feedback, while the latter is material that is still largely unpolluted 
(see metallicity profile below) and is condensing into gas clumps. 
Such gaseous clumps are found both at 50 and at 100kpc//i (co¬ 
moving) from the CM and represent counterparts of DLA systems. 
This is confirmed by the optical depth, rm, profile (third panel) 
which is boosted by high-density cold material. The optical depth 
markedly increases in correspondence to the center of mass of the 
halo (black solid line), reaching values of 10^. Additional peaks 
with thi > 10^ represent intervening cold material that gives 
rise to sub-DLAs or Lyman-limit systems. This is even more clear 
when looking at the off-center LOSs (red dashed and blue dotted 
lines) which are less impacted by the star formation episodes in 
the center of the halo and feature integrated optical depths up to 
thi ~ 10® —10^. We note here that si mpler HI modelings based on 
density cuts only ( such as the ones in iNagamine et al.ll2004[200^ : 
iTescari et al.ll200^ . reproduce similar trends for these off-center 
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Figure 14. Projected physical quantities extracted from three lines of sight close to the center of mass of the most massive halo at z = 6.594. Black solid line: 
LOS through the CM. Red dashed line: LOS at 50 comoving kpc/h from the CM. Blue dotted line: LOS at 100 comoving kpc/h from the CM. From top to 
bottom the panels show: total gas density (nu/cm^), gas temperature (K), H I optical depth (th i), H 2 density (nu/cm®), SFR (MQ/yr) and total metallicity 
(solar units). 


LOSs, although tend to overproduce rm values for particles near 
the central star forming region. 

In the fourth panel of Fig.[T4]we plot the H 2 density profile (result¬ 
ing from our non-equilibrium treatment) which is particularly inter¬ 
esting because of its enhancement corresponding to dense gaseous 
environments. Indeed, at the center of the halo the H 2 fraction 
experiences an exponential boost due to non-equilibrium catas¬ 
trophic runaway coolin j^ in dense collapsing sites and reaches 
peak values ~ 10 “^ with subsequent molecular-driven star for¬ 
mation episodes. Outside the halo’s zone of influence, molecular 
hydrogen does nol experience dramatic formation processes and 


We sU'ess that this is an out-of-equilibrium process taking place when 
densities increase and the gas passes from a loitering regime to an unstable 
regime without pressure support. 


hence the dotted and dashed lines never reach values as high as 
the black line. Nevertheless, local H 2 peaks can be easily identified 
with the dense material (see also first panel) at low temperatures 
(second panel) and with large rm (third panel) that has not yet un¬ 
dergone star formation episodes. 

This is clearly visible in the fifth panel where a spike in the SFR 
profile is present corresponding to dense H 2 -rich gas in the cen¬ 
ter of the halo, where stars are being formed actively (solid black 
line). Away from the central region (dashed and dotted lines), in¬ 
stead, SFR contribution is null - for the sake of clarity in this latter 
cases the SFR is arbitrarily set to logic SFR[M 0 /yr] = —3.2 in 
the panel (horizontal lines). 

The main effects of star formation are a global heating of the local 
environments and metal spreading into the surrounding medium. 
Stellar feedback is responsible both for H 2 destruction (see e.g. 
troughs of the H 2 profile in the fourth panel) and for the build-up 
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Figure 15. Simulated column density distribution functions, f{N), at 
2 = 2.995 (filled circles) and z = 4.992 (filled squares) and observa- 
tio nal determinations at z 3 from IProchaska et alj i2QQ5l) (diamonds) 
and iNoterdaeme etalJfeOlA (asterisks) with error bars. 


■ z = 4.992 
• z = 2.995 


^ Noterdaeme et al. (2012) 
^ Prochaska et al. (2005) 


^ ^ I 
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of metallicity gradients in the neighbourhood of star forming sites. 
In fact, the resulting metallicity profiles (last panel) are quite noisy 
and display significant amounts of heavy elements (up to ~ Zq) 
in correspondence with the broad temperature distribution around 
V « 1000 km/s shown in the second panel. Cold dense clumps far 
away from the central star forming regions are basically pristine. 
We will not go into the details of the individual elemental abun¬ 
dances, because they demonstrate shapes that are quite similar to Z, 
but we will perform a more statistical investigation. The discussion 
presented above has been based on three LOSs only and, although 
useful to shed light on the complex interplay between chemical and 
stellar feedback, it might suffer from low number statistics. 


4.2 Column density distribution function 


To have a statistically significant observable we compute the col¬ 
umn density distribution function, f{N), and make sure it is con¬ 
sistent with available observational determinations. This is a com¬ 
mon sanity check that we perform to be confident on how global 
properties of absorption systems are reproduced by our numerical 
calculations. 

We shoot 9000 random LOSs through the simulation box at each 
redshift and store the corresponding HI column densities. The col¬ 
umn density distribut ion function f{N) is then obtained, following 
previous works (e.g. lTvtleilll987tl . as the fraction of LOSs in the 
intervals (Wm, Nm + Nhi) ^ where A'hi denotes 

the HI column density and X the absoprption distance: 


X = 


Ho 

H{z) 


(1 -(- zfz. 


( 6 ) 


With these definitions, f{N)NX is the number of LOSs in the con¬ 
sidered X and HI column density bins. 

In Fig. [T^ we show f{N) extracted from our simulation box at 
redshift 2 = 2.995 (filled circles) and 2 = 4.992 (filled squares) 


with o bser vational determina t ions a t 2 ~ 3 from IProchaska et al.l 
( l2005h and iNoterdaeme et al.l ( l2012h . We find a generally decreas¬ 
ing behaviour consistent with available data for both DLAs and 
Lyman-limit systems and no evidence for strong redshift evolution 
in the epochs probed here. The drop at large Nm is rather steep, 
with a slope < —2. The trends for both 2 are essentially caused by 
the shorter cooling times at larger densities (see previous discussion 
in Sect. l3.2t . 

4.3 Statistics of thermal and chemical properties 

We extend the previous analysis to the cold gas clouds defined in 
Sect. 12.21 and studied in Sect. [3] For each object, we fix the origin 
of an orthogonal reference system on the CM, such that the 2 direc¬ 
tion points towards the observer. Then, for each object 17 LOSs are 
considered with displacements between 0 and ±15 kpc//i (physi¬ 
cal) from the center: 1 LOS along the 2 direction through the CM; 8 
LOSs radially displaced on the x axis (4 LOSs positively displaced 
and 4 LOSs negatively displaced); 8 LOSs radially displaced on 
the y axis (4 LOSs positively displaced and 4 LOSs negatively dis¬ 
placed). To randomize the viewing angle for each object a permu¬ 
tation of the X, y and 2 axes is performed. We take into account 
only absorbing systems with Nm ^ 2 x 10^° atoms cm~^ to have 
more directly comparable quantities with measurements. Since our 
implementation follows non-equilibrium species calculations and 
holds also in absence of pressure support it is possible to obtain 
reliable information on the gas thermal and chemical state and to 
explore the transition from Ft I atoms to H 2 molecules. 

In Fig. [T^ we show the mean temperature (top row), metallicity 
(middle row) and H 2 content (bottom row), all weighted by HI vol¬ 
ume density, as a function of HI column density. We consider the 
evolution with redshift from 2 = 6.594 (left column) to 2 = 4.992 
(middle column) and 2 = 2.995 (right column) for the resulting 
12587 LOSs (680 at 2 = 6.594, 2907 at 2 = 4.992 and 9000 at 
2 = 2.995). In each panel of the figure the solid orange line shows 
the median values. The blue triangular points refer to LOSs passing 
through the CM, while the black squares refer to off-center LOSs. 
At a fixed HI column density we notice some scatter in the temper¬ 
atures (especially at low Nm), even if the median remains constant 
around 10^ K across the whole column density range. This scat¬ 
ter reflects the different environmental conditions of the absorb¬ 
ing systems and the effects of intervening hot material between the 
observer and the targeted structures. Moving from 2 = 6.594 to 
2 = 2.995, the majority of LOSs report values of the H I-weighted 
median temperature below ~ 3 x 10^ K (i.e. 4.5 dex), consistently 
with our discussion in Sect.[3 The rare LOSs at logic [T/K] > 4.5 
(1 at 2 = 6.594 and 8 at 2 = 2.995) intercept parcels of gas just 
experiencing bursts of star formation. The LOSs with lower temper¬ 
atures are instead passing through dense gas dominated by molec¬ 
ular (H 2 and HD) cooling and undergoing run-away collapse. 

In the second row of Fig. [T^ we present the redshift evolution of 
the H I-weighted total metallicity (in solar units) along the different 
LOSs. At a fixed Nm the metallicity value is severely scattered, 
with a small fraction (of the order of ~ 10 per cent) of values at 
Z < 1O-^Z0. We note that, given these numbers and the fact 
that only a bunch of metal-poor DLAs are currently known, it is 
statistically difficult to draw conclusions on such low-metallicity 
objects and further observational data are definitely needed. It is 
also difficult to measure metallicities at redshift 2 ~ 6.6, because 
of the lack of informatio n about the HI column density due to the 
I Gunn & PetersonI d 19651) effect. The median metallicity is roughly 
constant at Z < 10 “^Zq and increases by one order of magni- 
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Figure 16. Physical quantities of simulated DLAs. Gas mean temperature (top row), total metallicity (middle row) and H 2 mean density (bottom row) as a 
function of HI column density. We consider three different redshifts: z = 6.594 (left column), 2 : = 4.992 (middle column) and 2 = 2.995 (right column). 
In each panel, the solid orange line shows the median of the corresponding physical quantity. Blue triangles refer to LOSs passing through the center of mass, 
while black squares refer to off-center LOSs. 



Figure 17. H I-weigthed metallicity distributions at redshift 2 = 6.594 (left panel), 2 = 4.992 (central panel) and 2 = 2.995 (right panel). 


tude (> 10~^Zq) in high-density regions due to ongoing metal 
pollution between 2 = 6.594 and 2 = 2.995. Already at earlier 
epochs (2 = 6.594), the median Z shows a large metallicity spread 
as a result of the initial stages of metal pollution: the first phases 
of this process in the denser cores of primordial structures cause 
very inhomogeneous enrichment patterns. At 2 = 4.992 the LOS 
metallicities remain around 10“^ Zq, despite some outliers below 
10-4Zq. 

At 2 = 2.995 the metallicity distribution of the LOSs peaks still at 
~ 10“^ Zq, while Z < IQ-"* Zq data account for ~ 10 per cent of 


the cases. Interestingly, we observe a steep increase in the high-Nn i 
regions caused by ongoing pollution from the numerous low-mass 
stars that can now die as SN la (see discussion in Sect. 13.3.^ . 

For sake of completness, in Fig.[T3we show the corresponding HI- 
weighted metal distribution functions at 2 = 6.594, 2 = 4.992 and 
2 = 2.995. As previously mentioned, there is little redshift evo¬ 
lution, although at earlier times Z values are more sharply peaked 
around 0.01 Zq (2 = 6.594) and later on between 0.01 — 0.1 Zq. 
Finally, in the bottom row of Fig. [H we plot the FI 2 content as a 
function of HI column density. Larger-Nm structures sample more 
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gas and hence have more chance to host H 2 formation, however, 
it is not simple to state whether H 2 correlates with Nhi along the 
LOS. The trends feature evident temporal and local variations of 
H 2 abundances accompanied by a broad spread. The denser regions 
host the first star formation events and metal spreading episodes. 
Star formation feedback is then responsible for molecule dissoci¬ 
ation in the surrounding environment. The resulting behaviour is 
therefore a consequence of the increasing star formation activity 
with cosmic time. In fact, its feedback dissociate molecules and af¬ 
fect the slope of the HI-H 2 relation. 

As a final remark, we stress that both shielding corrections and 
weighting scheme could have impacts on the final results. We have 
checked that w i th the altern ative shielding prescript ions by e.g. 
iNagamine et al.l ( |2004[200% and iTescari et alj l l2009l) we get typ¬ 
ically lower scatter in the HI-H 2 relation, due to t he lack of tem- 

S ieratu re effects, which are instead accounted for bv iRahmati et al.l 
2013ll . Different weigh ting schemes can i ntroduce further system- 
atics. Ongoing studies jRubin et alj|2014h or future observational 
campaigns might be able to test LOS modeling. 

Thermal properties for off-center LOSs (black squares) and for 
LOSs passing through the CM (blue triangles) do not differ signif¬ 
icantly at any redshift. By looking at the latter case we can confirm 
that, in the central regions, high H 2 values drive non-equilibrium 
runaway gas collapse and star formation. Subsequent feedback 
causes gas evacuation, metal spreading and molecule dissociation 
with a remarkable decrement in the H 2 fractions. Very recent ob- 
servational analyses of 2 2 — 3 data bv iNoterdaeme et aP l l2015tl 

seem to agree with the expected large scatter in H 2 abundances. 


5 CONCLUSIONS 


We have employed cosmological N-body hydrodynamical chem¬ 
istry calculations following gas atomic and molecular collapse and 
heavy-element production from stars with different masses and 
metallicities during the first few Gyr of the Universe, in the red- 
shift range 2 < 2 < 7. Our detailed non-equilibrium chemistry 
integration of e“, H, H+, H~, He, He+, He++, D, D+, H 2 , H+, 
HD, HeH^ abundances and stellar evolution treatment with indi¬ 
vidual yields for different species (He, C, N, O, Ne, Mg, Si, S, Ca, 
Fe) allowed us to investigate the chemical patterns of intermediate- 
redshift galaxies and compare the predictions of our calculations to 
observational measurements. 

Throughout this work we have assumed a ACDM scenario. We 
warn the rea der, though, that in s ome particular cases (e.g. warm 
dark matter, iMaio & Viell l2015l) the background cosmological 
model can influence the baryonic-structure evolution, mostly at 
higher redshift. This can play a role during the very beginning of the 
onset of star formation, but we checked that typically star formation 
and feedback mechanisms tend to dominate and alleviate possible 
discrepancies in the epochs of interest here (see Maio et ahl 200d; 
iMaio & lannii^ I 2 OI II: IPace & Maid I20l4 iMaio & Viell l2015h . 
Similarly, the exact numerical parameters adopted for the initial 
mass functions, wind prescriptions, initial-condition gas velocities 
and related issues are likely to ind uce some changes , which are ex¬ 
pected to be minor (see more in iM aio et al ]|201 ll : ICamnisi et al.l 
I 2 OIII: iTescari et alj|2014l : 1^ et al.ll2015h . Uncertainties might de¬ 
rive from the lack of a definitive treatment of diffusion mechanisms, 
which, despite several attempts, still remain an unsolved problem 
in astrophysics. 

Theoretical yields are affected by a plethora of uncertain physi¬ 
cal processes in stars (such as explosion mechanisms, differential 


rotation, the initial composition, magnetic fields, nuclear reaction 
rates, etc.) all of which can influence the final ratios. Neverthe¬ 
less, we find agreement between our theoretical expectations and 
observational determinations. We note that the detailed values of 
metal yields for individual elements are not supposed to change 
gas hydrodynamics significantly, because metal content and cool¬ 
ing in the regimes explored here are usually dominated by oxygen, 
for which there are well determined metallicity-dependent yields 
(e.e. lFrancois et al.ll2004) . 

Physical processes in the inter-stellar medium are modelled by 
means of subgrid effective models, w hich have also be en tested 
through simulations of isolated objects jMaio et al.ll2013h . The ef¬ 
fects of changing particular model parameters have been shown not 
to be dramatic for star formation and enrichment processes, as long 
as gas cooling capabilities (which have been treated in a quite de¬ 
tailed fashion in this work) are properly addressed for both atomic 
and molecular phases. 

We do not expect statistically significant changes when consider¬ 
ing different IMFs at early times due to e.g. pristine (popIII) gen¬ 
erations. Indeed, previous works llTornatore et al.l2007f:lMaio et al.l 
20ld:IWise et aljboiTl : [Muratov et al.ll20I3l:|Pallottini et al.ll20l4 


Hirano et alj|2015[) have shown that these should be cosmologi- 


cally subdominant at any redshift, almost irrespectively of their 
mass spectrum. Isolated regions with pristine pockets of popIII star 
forming gas can survive down to very low redshift, however, they 
are not supposed to contribute significantly to cosmic SFR densi¬ 
ties (Fig. IS, neither in the cas e of massive pop III stars nor in the 
case of low-mass popIII stars iMaio et al.ll201 ll) . Further more, re¬ 
cent a nalyses based on measured GRB host metallicities iMa et al.l 
[ 2 oTi) seem to exclude massive popIII star formation episodes for 
current detections at 2 > 4.7. 

We find that redshift ~ 2 — 7 galaxies have baryonic masses up to 
10^^ M 0 with typical stellar masses as low as 10® M©. Struc¬ 
tures at intermediate redshifts have typical SFRs around 10~^ 
and ~ 20MQ/yr and a very broad SSFR range of ~ 10“^ — 
10^ Gyr“^ with typical values of ~ 10Gyr“^ at 2 > 6.5 and 
more quiescent ones around ~ 0.1 — 1 Gyr“^ at 2 < 3 (Fig|5](. 
The combined effects of feedback and environment can cause large 
excursions in the stellar fraction with typical values extending from 


SFRs (< 10“^ M 0 /yr) and stellar fractions (< 10“^) and rep¬ 
resent objects that are just b uilding up their ste llar mass (consis¬ 
tent with observational hints: IFvnbo et al.ll2008h . Our results sug¬ 
gest that metal-poor DLAs can be transient objects formed by con¬ 
densation of warm material (~ 10“* K) and captured in the evo¬ 
lutionary stage when cooling processes are lowering the kinetic 
temperature, allowing the gas to condense to higher overdensi¬ 
ties (Figs. [6l7t . This explains why, despite molecular (H 2 ) frac¬ 


tions might be locally large, their avera ge values are usua l ly ver y 
scattered, cons i stently with data by e.g.lN oterdae me et al.l ( l2008h : 
ISrianand et alj feoioh : lAlbomoz Vasguez et alj l l2014h . To explic¬ 
itly demonstrate this evolutionary phase of metal-poor DLAs, we 
have traced a few simulated gaseous clumps over the redshift inter¬ 
val 2 ~ 2.995 — 6.594 and have shown their large-scale contrac¬ 
tion. On the other hand, we have additionally checked the possi¬ 
ble fate of metal-poor DLA candidates formed at high redshift and 
found that they are likely to be partially or totally destroyed due to 
feedback from in situ star formation or tidal interactions. 

In the comparisons with available data, the metallicities and molec¬ 
ular fractions derived from our numerical calculations are in agree¬ 
ment with the observations. 

By tracking the cosmic chemical evolution for some of the most 
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common elements and by testing the resulting abundance ratios 
(such as [Fe/H], [O/H], [C/0], [Si/0], [0/Fe], [Si/Fe], [N/0]) 
against observational measurements (Figs. mm, we find that aver¬ 
age trends and dispersion of simulated systems agree with observed 
systems from z ~ 2.995 to z ~ 6.594. The ratios we consider 
do not present strong redshift evolution, since metal enrichment is 
dominated by massive stars at these epochs with some contribu¬ 
tion from intermediate-mass stars and SN la. In all cases, mean 
values for e.g. [C/0], [O/F e], [Si/Fe] behave within the error bar s 
of observational estimates dBecker et al.ll20l3 : ICooke et alj|2015h . 
The late-time Fe enhancement from long-lived SNe la decreases 
[O/Fe] and [a/Fe] ratios consistently with stellar evolution time- 
scales (Figs. imn). 

An interesting feature of our analysis is the existence of consider¬ 
able scatter along the theoretical trends over the redshift intervals 
we probe: typical standard deviations vary from a fraction of one 
dex up to more than half of a dex. We interpret these spreads as due 
to simultaneous enrichment by stars at different stellar stages or 
belonging to different stellar populations. This is a remarkable de¬ 
viation from the simple-stellar-population scenario and highlights 
the limits of homogeneous models. 

From an observational point of view, the effects of SN la going off 
at z < 4 — 5 should appear in the broadening of measured line 
widths at later times as a consequence of the establishment of a 
more and more turbulent regime. 

Simulated spectral observations (Fig. IS} allow us to stress the con¬ 
nection between gas, H 2 , star formation and the complex interplay 
between chemical and stellar feedback in the presence of cold opti¬ 
cally thick molecular-rich clumps. Our synthetic spectra statistics at 
z-2-7(Fig.[l6j suggest that the median intervening gas temper¬ 
ature evolves little during structure growth after reionization. The 
increase in metallicity from high to lower z and the large scatter 
observed in H 2 abundances are mainly due to the increasing star 
formation activity (Figs. [2l5) in increasingly larger objects, after 
the primordial bursty phases. 

These findings demonstrate the importance of baryon chemistry to 
understand the properties of cosmic structures at early and interme¬ 
diate epochs and its power to unveil the origins (SN II, AGB and 
SN la) of heavy elements and of their trends during cosmological 
evolution. 

Despite disentangling theoretical expectations and observational 
features might be challenging, enrichment of metal-poor gas could 
hide indirect signatures of previous stellar generations. Future stud¬ 
ies of elemental ratios as expected from different stellar populations 
within 3-dimensional numerical simulations will shed light on the 
properties of low-Z star formation and on the existence of primor¬ 
dial popIII stars. 
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APPENDIX A: DUST ESTIMATES 

The sudden appearance of large amounts of dust in early epochs 
can be explained b y cooling and condensation of heavy elements 
Jpeder et al.l Il966h in massive core-collapse SN ejecta. These, 
thanks to the short SN lifetimes and the high melting tempera¬ 
tures of the principal refractory elements (Si, Mg, Fe, Ca, Ti, Al, 
etc.), act as major sources of dust production th at can be addi¬ 
tionally reprocessed in th e inter-stellar m edium jYozin & Bekkil 
|2014 iMancini et alj|201^ : lMattsson|l2015h . Investigations of ex¬ 
ternal galaxies and Magellanic Clouds have suggested a tight de¬ 
pendence of the dust-to-gas ratio, D, on the local metallicity, Z. A 
stringe nt correlation ha s been obtained via sub-millimeter analyses 
by e.g. iGalametz et al.l (1201 ih an d has been suppor ted by analyti¬ 
cal models of dust grain evolution dFeder et al.l 19661) . The resulting 
trend shows an increasing D for increasing Z, D oc Z^'^, over a 
broad range of Z. Although averaged over the many complex and 
still poorly understood processes that rule formation and destruc¬ 
tion of dust grains (e.g. the role of different stellar sources, metal 
yields, local star formation, radiative transfer, sticking coefficients, 
grain size distribution, grain destruction/growth timescales in the 
inter-stellar medium, reverse-shock effects, underlying models for 
grain evolution, etc.), the D{Z) behaviour can be very precious to 
get hints on the typical fraction of metals condensed onto dust in 
different metallicity regimes. Environmental effects can be surely 
importan t, as also testified by the large error bars for the data re¬ 
gression dGalametz et al.llToilh . but the general trend is fairly es¬ 
tablished. Thus, the empirical D{Z) relation represents a cheap and 
fruitful tool when applied to simulated cosmic structures, since it 
implicitly includes and constrai ns all the relevant unkn owns. Here, 
we apply the relation fitted by IGalametz et al.l d201lh to infer the 
dust mass implied by the metal content of each object in our data 
sample. 

In Eigs. l A1 l and lA2l we show the expected dust mass as a function of 
gas mass, average molecular fraction, stellar mass and specific SER 
for the whole data sample at redshift z = 6.594 and z = 2.995, 
respectively. 

We observe already at higher z a well established increasing trend 
of dust mass, Ma, with gas mass. The amounts of dust can be as 
large as ~ 10^ — 10® Mq already in the first billion years and the 
corresponding average D value can reach ~ 0.007. The effects of 
metal spreading is reflected in the broad Ma distribution which de¬ 
parts from the average behaviour (empty diamonds), mostly in the 
low-mass end. When looking at the trend with molecular content 
we do not see any particular correlation, due to the dominant role 

http://nf.nci.org.au 

https://edward-web.hpc.unimelb.edu.au/users/ 
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Figure Al. Dust mass estimates as function of gas mass (top left), average 
molecular fraction (top right), stellar mass (bottom left) and specific SFR 
(bottom right) at 2 = 6.594. The average dust-to-gas ratio, D, is quoted in 
the top left panel, with corresponding expectations (red empty diamond). 


Figure A2. Dust mass estimates as function of gas mass (top left), average 
molecular fraction (top right), stellar mass (bottom left) and specific SFR 
(bottom right) at 2 = 2.995. The average dust-to-gas ratio, D, is quoted in 
the top left panel, with corresponding expectations (red empty diamonds). 


of metallicity. The behaviour with stellar mass, is instead, more 
clear with peak values reached for M* ~ 10® — 10^® Mq. Interest- 
ingl y, this is consistent with recent observations at z ^ 6.5 — 7.5 
(e.g. lSchaerer et alj2015l : IWatson et al.l2015l : lMaiolino et al.l20T^ 

whose data suggest normal galaxies with stellar masses M* ~ 
10® — 10^® Mq and dust masses Ma 10^ — 10® Mq. In terms 
of SSFR, more quiescent objects host typically larger dust content, 
as a hint of their previous activity. 

At later times (z = 2.995) the picture shown in Fig. IA2l is simi¬ 
lar, although the typical dust-to-gas ratio increases by a factor of 2, 
to D 0.014, and the dust fraction is comparable to metallicities 
(see main text). 

From a qualitative point of view the over-linear dependence of 
D (X ® helps understand why dust issues should be more rele¬ 
vant in more enriched gas. Finally, we note that, despite the crude 
approximations, our approach based on the dust-metallicity relation 
performs relatively well when compared to data or to more sophis¬ 
ticated models that have to include unknown free parameters. 
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